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In the Minimal Supersymmetric Model (MSSM) the hitherto unconstrained flavor mixing between 
top-squark and charm-squark will induce the flavor-changing neutral-current (FCNC) interaction 
between top quark and charm quark, which then lead to various processes at the next generation 
linear collider (NLC), i.e., the top-charm associated productions via e + e~, e~7 and 77 collisions as 
well as the top quark rare decays t — * cV (V — g, 7 or Z). All these processes involve the same 
part of the parameter space of the MSSM. Through a comparative analysis for all these processes 
at the NLC, we found that the best channel to probe such SUSY-induced top quark FCNC is the 
top-charm associated production in 77 collision, which occurs at a much higher rate than e + e~ or 
e~7 collision and may reach the detectable level for some part of the parameter space. Since the 
rates predicted by the Standard Model are far below the detectable level, the observation of such 
FCNC events would be a robust evidence of new physics. 
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I. INTRODUCTION 

The study of flavor-changing neutral-current (FCNC) interactions has been playing an important role in testing 
the Standard Model (SM) and probing new physics beyond the SM. As the most massive fermion in the SM, the top 
quark is naturally regarded to be more sensitive to new physics than other fermions. fn the SM, due to the GIM 
mechanism, the top quark FCNC interactions are absent at tree level and extremely small at loop levels. In some 
new physics models beyond the SM the top quark FCNC may be significantly enhanced. Searching for the top quark 
FCNC would be a good probe for new physics. 

Intensive activities to explore the top quark FCNC couplings have been undertaken in recent years. On the 
experimental side, the CDF and DO collaborations have reported interesting bounds on the FCNC top quark decays 
[1] from Run 1 experiment and will tighten the bounds from the on-going Run 2 experiment. On the theoretical 
side, various FCNC top quark decays and top-charm associated productions at high energy colliders were extensively 
studied in the SM [2,3], the Minimal Supersymmetric Model (MSSM) [4-7] and other new physics models [8]. These 
studies showed that the SM predictions for such top quark FCNC processes are far below the detectable level and the 
MSSM can enhance them by several orders to make them potentially accessible at future collider experiments [9]. 

Due to its rather clean environment, the next generation linear collider (NLC) will be an ideal machine to probe 
new physics. In such a collider, in addition to e + e~ collision, we can also realize 77 collision and e~7 collision with the 
photon beams generated by the backward Compton scattering of incident electron- and laser-beams [10]. The SUSY 
induced top-charm FCNC will give rise to various processes at the NLC, i.e., the top-charm associated productions 
via e + e~, e~j and 77 collisions as well as the top quark rare decays t — » cV (V = g, 7 or Z). It is noticable that 
some of these processes, like the top-charm associated productions in 77 or e~7 collision, have not been studied in 
the framework of the MSSM. The production in 77 collision may be more important than in e + e~ collision studied 
in the literature [7]. The reason is twofold. Firstly, the process 77 — > tc is a good probe of new physics because it 
is essentially free of any SM irreducible background [11]. Secondly, unlike the process e + e~ — > tc, which is s-channel 
suppressed in high energy collisions, there are t- and u-channel contributions to 77 — > tc and thus its cross section 
may be much larger at the NLC. It is also noticable that all these FCNC processes at the NLC involve the same 
part of the parameter space of the MSSM 1 . Therefore, it is necessary to perform a comparative analysis for all these 
processes to find out which process is best to probe the top quark FCNC. This is the aim of this article. 

This paper is organized as follows. In Sec. II, we discuss the mixing between top-squark and charm-squark and 
derive the FCNC interaction Lagrangian in SUSY-QCD. In Sec. Ill, we calculate the tc productions in 77 and e~7 



1 Since the FCNC decay t — > ch involves more parameters, we do not include it in our analysis. 



1 



colllisions induced by such FCNC SUSY-QCD interactions. Numerical results for these productions at the NLC are 
given in Sec. IV, with the comparison to tc production in e + e~ collision and various FCNC top quark decays. Finally, 
a brief discussion and conclusion can be found in Sec. V. 



II. FLAVOR MIXING OF TOP-SQUARK AND CHARM-SQUARK AND THE INDUCED FCNC 

Many popular SUSY models predict the flavor mixings of sfcrmions. For the squark sector, despite of the possible 
strong constraints on the down type squark flavor mixings from the low-energy experimental data 2 , the mixings 
between top-squark and charm-squark are subject to no strong low-energy constraints [13]. 

Such a mixing between top-squark and charm-squark is well motivated in low-energy supergravity models (SUGRA) 
[14]. In these models squark mass matrices are diagonalized simultaneously as with the quark matrices at the Planck 
scale. But when the terms evolve down to low energy, this diagonality is violated by radiative corrections. Due to 
large top quark mass, at low energy it is found [15] that the mixing between top-squark and charm-squark may be 
significant. Note that only cl mixes with top-squark while cr does not in the approximation of neglecting charm 
quark mass. As shown in [15], the mixing between cl and II is most likely to be large, which is proportional to a 
sum of some soft masses. 

Motivated by the above arguments, in our analysis we assume the existence of the mixing between cl and <l and 
parameterize the mixing as 5lMqMq i , where Mq (Mq i ) is the soft mass parameter for left-handed squark of third 
(second) generation and 8l is a dimensionless parameter representing the mixing strength. 6l can be calculated in 
terms of other parameters in a given model like mSUGRA [15]. But in our calculation we retain it as a free parameter 
in the range of ~ 1 [13]. 

Considering the mixing between II and Ir and neglecting the mixing between cl and cr, we obtain the 3x3 
mass-square matrix in the basis (II, tR, £l): 



( t* L t* R c* L ) 



( M^+m 2 + D L m t X 

m t X M? • n,j ■ D R | | t R | (1) 



V S l MqM^ 




where Dl = m 2 z cos(2/3)(| — §sfy), Dr = cos{2j3)s 2 v , X = A t + /Ltcot/3 and M?. is the soft-breaking mass term 

for right-handed top-squark. A t is the coefficient of the trilinear term H2QU in soft-breaking terms, fi is the mixing 
mass parameter between Hi and H 2 in the superpotential and tan j3 — v 2 /v\ is ratio of the vacuum expectation values 
of the two Higgs doublets. 

The mass-square matrix in Eq.(l) can be diagonalized by a 3 x 3 unitary matrix V , which rotates the interaction 
eigenstates (II, tR, Cl) into mass eigenstates 91,2,3- Such a rotation results in the FCNC in both the weak interaction 
sector and the strong SUSY-QCD sector. Of course, the latter will be dominant. So in our analysis we only consider 
the FCNC in SUSY-QCD sector, given by 

Cfcnc = -V2g s T a [§ a (V&P L - V? 2 P R )tq* + ], a V^P L cq*\ , (2) 

where P L . R = (1 T 7s)/2. 



III. CALCULATIONS 



The SUSY-QCD FCNC in Eq.(2) will induce the top-charm associated production in 77 collision, as shown in 
Fig. 1. Throughout our calculations the charge conjugate production channel, i.e., the production of tc, has also been 
included. Neglecting the charm quark mass, the amplitude of this process is given by 

M - 2iaa s Q 2 c C F u t Y^P L v c e„(Ai)e„(A 2 ) (3) 

with Qc = 2/3, C F = 4/3 and 



2 The down type squark flavor mixings in large tan/? limit could enhance the FCNC B decays by several orders [12] and will 
face tests in the on-going B-factory experiments. 
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c 9 g^ + cioW + diPtPth + ci2P<tPcfa + dzPtPc^ + diPtPch + c ^Pt K 



2/3- 



(4) 



Here ki, 2 denote the momentum of incoming photons and p t , c the momentum of outgoing top and charm quarks. The 
coefficients Cj can be obtained by a straightforward calculation of the diagrams shown in Fig. 1 and their expressions 
are presented in the Appendix. 
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FIG. 1. Feynman diagrams contributing to the process 77 — > tc in SUSY-QCD. 



We checked that our results satisfy the Ward identity, k^T^ = and k^T^ = 0. We also checked that all ultraviolet 
divergences canceled in our results due to the unitary of V, which is essentially guaranteed by the renormalizability 
of the MSSM. Note that when t or u approaches 0, which is physically permitted, Cj tend to be very large. Indeed, 
this is the advantage of such a process over e + e~ — > tc which has only s-channel contribution and thus is suppressed 
at high energy colliders. 



3 



The tc production in e~7 collision proceeds through the process e~7 — ► e _ 7*7 — » e~tc, where the 7-beam is 
generated by the backward Compton scattering of incident electron- and laser-beam and the 7* is radiated from e~ 
beam. The subprocess 7*7 — ► tc has the same Feynman diagrams shown in Fig. f. In our calculation we use the 
Weizscicker- Williams approximation [16] which treats 7* from e~ beam as a real photon. Thus the process can be 
approximated by the simpler fusion reaction 77 — > tc and its cross section is given by 

&e~y^>etc{Se-y) / dx Pj/ e (x : E e ) (7 77 ^£c( s 77 — xs ej): (5) 

J (m t +m ) 2 /s e7 

where P 1 / e {x, E e ) is the probability of finding a photon with a fraction x of energy E e in an ultrarelativistic electron 
and is given by [16] 



a /l + (l-ar) 2 /, E e 1\ x / ,2 ^ \ (2 - a;) 2 , /2-2ar M 



Note that there are also intermediate Z-boson contribution for the process — > etc. However, such contributions 
are suppressed by the probability function of finding a Z-boson in an ultrarelativistic electron [17] and can be safely 
neglected. 

For both 77 collider and collider, the photon beams are generated by the backward Compton scattering of 
incident electron- and laser-beams just before the interaction point. The events number is obtained by convoluting 
the cross section with the photon beam luminosity distribution. For 77 collider the events number is obtained by 



^77— >tc 



/ ^\/^77 1 , ^"77—^*0(^77) — £ee ^"77— >tc (See ) > (7) 

J V ^77 



where d£ 77 /d v /s^ 7 " is the photon beam luminosity distribution and cr 77 _^i 5 (s e e), with s ee being the energy-square of 
e + e~ collison, is defined as the effective cross section of 77 — ► tc. In optimum case, it can be written as [18] 

rx malc rx max ^ z 2 

(T 77 _ te (s ee ) = / _2zdz <T 77 ^ tE (s 77 = z 2 s ee ) / — F l/e {x) F l/e {— ), (8) 

J (m t +m c )/ ^/s cc Jz 2 /x maus x x 

where F 7 / e denotes the energy spectrum of the back-scattered photon for unpolarized initial electrons and laser photon 
beams given by 

/ x 1 (■> 1 4x 4x 2 \ 

The definitions of parameters £, D(£) and x max can be found in [18]. In our numerical calculation, we choose £ = 4.8, 
£>(£) = 1.83 and x max = 0.83. 

For the collider the effective cross section of e-f — > etc is defined as 



1 r , d£ £ 

5e7 



<7e7— >etc(^ee) ~ / d^/ S e7 <J e ^—* e tc { s e~f ) 



= _2zdz ir 77 _ >tE (s 77 = 2; 2 s ee ) / — P^/ e {x,E e )F l/e {—). (10) 

J (m t +m c )/^/s^ Jz 2 /x ma!c x x 

The process e + e~ — > tc and the top quark rare decays t — > cU (V = 7, Z, g) were already calculated in the literature. 
Here for comparison we recalculated all of them. The lengthy expressions are not presented here. 



IV. NUMERICAL RESULTS 



The relevant SUSY parameters all the top-charm productions and top quark rare decays are Sl, X (= A t + /x cot /3), 
Mg 7 Mq, Mjj, Mq X and tan/3. To find out typical magnitudes of these processes, we scan over these SUSY parameters 
by requiring 5 < tan/3 < 50, m, > 86.4 GeV [19], m g > 190 GeV [19], rriQ 1 ~ 1 TeV and restricting other soft mass 
parameters to be of sub- TeV scale. 

Our scan results are plotted in Figs. 2 and 3. For the productions one sees that the typical values are (1(77 — > 
tc) ~ 10~ 2 fb and cr(e7 — > etc) ~ 10~ 3 fb and in optimum cases (7(77 — > tc) ~ 0.7 fb and cr(e7 — > etc) ~ 0.04 fb. 
While for the process e + e~ — > tc, the cross section can only reach 0.02 fb, one order of magnitude lower than that of 
77 — > tc and comparable with that of e7 — > etc. The reason for this, as we pointed out before, is due to the t-channel 
enhancement for 77 — ► tc and s-channel suppression for e^e - — > tc. For the top quark rare decays we see that the 
optimum values are Br(t — ► eg) <~ 10 -5 , £?r(t — » cZ) ~ 10 -7 and _Br(t — » C7) <~ 10~ 7 , which agree with previous 
results [5]. 
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FIG. 2. The scattered plot of SUSY-QCD contribution to processes 77 — > te, e7 — > etc and e e — ► tc for ^/s ee = 500 GeV. 




FIG. 3. Same as Fig. 2, but for FCNC top quark rare decays. 



Comparing with predictions in other models, such as R-parity violating model [20] and the type III two Higgs 
doublet model [21], we find that the optimum value of cr(77 — ► tc) in SUSY-QCD is larger. Note that our results for 
a(e + e~ — > tc) are quite different from those in Ref. [7] where the optimum value of er(e + e~ — > tc) can reach 0.1 fb. 
The reason for the difference is that in Ref. [7] a non-unitary flavor changing matrix is used. During our scan, we 
find that, only in case of a light top-squark with mass of 100 ~ 200 GeV and a large mass splitting of a few hundreds 
GeV in squark spectrum can relative large cross sections of these processes be obtained. This is due to a GIM-likc 
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cancellation which occurs in the limit of degenerate squark masses as a result of the unitarity of the rotation matrix 
V. 

The SM and MSSM predictions for the top-charm FCNC processes are summarized in Table I. Note that since 
the SM prediction for the process 77 — > tc has not been done in the literature, the SM value shown in the table is 
calculated by us. This is an arduous work. There are 66 diagrams contributing 77 — -> tc in gauge and often the 
contribution from a single box diagram can be decomposed into dozens of terms. In order to avoid artificial error, 
we use FormCalc [22] in our calculations. From Table I, we sec that the MSSM can enhance the cross section of 
top-charm associated production processes by several orders. 



Table 1: Theoretical predictions for the top-quark FCNC processes. SUSY-QCD predictions are the maximum values. The 
collider energy is 500 GeV for productions. 





SM 


SUSY QCD 


(7(77 — > tc) 


0(l(T a ) fb 


©(lCr 1 ) fb 


er(e~7 — * e~tc) 


0(l(T y ) fb 


O(10- a ) fb 


<r(e + e~ — > tc) 


e>(icr iU ) fb 


ollO-' 4 ) fb 


Br(t -> eg) 


©(lCr 11 ) 


O{10' b ) 


Br(t -> cZ) 


olio- 1 -') 


O(10-') 


Br(t — > cry) 


0(icr ia ) 


O{10-') 



The behaviours of the SM cross sections versus the collider energy are shown in Fig. 4 for three different production 
processes. We see that the cross section of e + e~ — > tc drops quickly with the increase of collider energy due to the 
s-channel suppression. 
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FIG. 4. Cross sections of top-charm associated production processes versus collider energy E c 



/s ee in the SM. 



Let us focus on the most important process 77 — > tc and study some of its features. For example we fix tan/3 = 30 
and Mq 1 — 1000 GeV, and choose Mq = 600 GeV, Myy = 400 GeV for scenario I (relatively light squarks for the 
third- family) and Mq = 800 GeV, Mq = 600 GeV for scenario II (relatively heavy squarks for the third- family) . 
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X = 800GeV 
E cm = 500GeV 
m § = 200GeV 



FIG. 5. The SUSY-QCD correction to o" 77 ^js as a function of 5l- The solid curve corresponds to scenario I (relatively light 
squarks for the third-family) and the dashed curve corresponds to scenario II (relatively heavy squarks for the third-family). 
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FIG. 6. Same as Fig 5, but for (7(77 — > tc) as a function of X (= At + ficot /3). 



The cross section of 77 — > tc as a function of 5^ is shown in Fig. 5. This figure shows that the cross section is 
enhanced dramatically with the increase of Sl- Clearly this is due to the fact that large 5l will not only enhance the 
mixing between top-squark and charm-squark but also enlarge the mass splitting between the squarks. In Fig. 6 we 
show the dependence of (7(77 — * tc) on X. We see that (1(77 — > tc) increases as X gets large. The reason is that 
large X enhances the mass splitting between the top-squarks and thus lead to a weak cancellation between different 
Feynman diagrams. 
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FIG. 7. Same as Fig 5, but for the dependence on luino mass. 




The dependence of 17(77 — * ^c) on m,g and y/s is shown in Fig. 7 and Fig. 8, respectively. Fig. 7 shows that the 
cross section drops when gluino becomes heavier, showing the decoupling effect of the MSSM. 

From Figs. 5-8 one finds that the cross section in scenario I is always larger than that in scenario II. The reason 
is that scenario I has the relative light squarks for the third-family which lead to a larger mass splitting in squark 
spectrum. 
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V. SUMMARY AND CONCLUSION 



We have performed a comparative analysis for the SUSY-induced top-charm FCNC processes at the NLC. These 
processes include the top-charm associated productions via e + e~, e~j and 77 collisions as well as the top quark rare 
decays t —* cV (V = g, 7 or Z), all of which involve the same SUSY parameters. 

In the production channels the tc production in 77 collision was found to occur at a much higher rate than e+e~ 
or e~7 collision. In some part of parameter space, the production rate of 77 — ► tc can reach 0.7 fb. This means 
we may have 70 events each year for the designed luminosity of 100 fb _1 /year at the NLC. Since the SM value of 
the production rate is completely negligible, the observation of such tc events would be a robust indirect evidence of 
SUSY. 

Note that in practical experimental searches of such tc events, a careful study of backgrounds is needed. To 
efficiently suppress the backgrounds, the signal events may be hurt by 50% or more [11,23]. Using polarization at a 
linear collider can help to suppress the background, as analyzed in the second reference of [9]. 

In the rare decay channels the t — > eg was found to have the largest branching ratio which can reach 10~ 5 , in 
agreement with previous studies. Although 10 5 tt events could be produced at the NLC each year, studies [24] showed 
that the sensitivity to such rare decays can only reach 5 x 10~ 4 . So the MSSM prediction for the branching ratio of 
10~ 5 is too low to be accessible at the NLC unless the designed luminosity can be further upgraded. 

Therefore, we conclude that, to probe the SUSY-induced FCNC top quark interactions at the NLC, tc production 
in 77 collision is the best channel. 

Note that in our analysis we assumed flavor mixing occurs between cl and II, which is favored in SUGRA models. 
Other kinds of mixings like the mixing between cr and t^ [25,26] may also be of phenomenological interest. However, 
we believe our conclusion will be quite model-independent since our results reflect the basic features of the processes. 
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APPENDIX 



Before presenting the explicit form of , we define the following abbreviations 

8= {kx + k 2 f = (pt+Pcf, (11) 
t = (pt-k 2 ) 2 = (fci-;p c ) 2 , (12) 
u= {pt-h) 2 = {k 2 - Pc f, (13) 
3 

V 13 B? = £ V^Bii-pum^ms), (14) 
0=1 
3 



V 13 B\ = V 10 V^B t ( Pc , m 0p , m- g ), (15) 
0=1 
3 

VuBi = E VipV^Bii-pt + k 2 ,m O0 ,m~ g ), (16) 
0=1 
3 

VisCfj = V^V^djikuPc - ki,m 0p ,m &p ,m 5 ), (17) 
0=1 
3 

VwCfj = Y ^l/jVjgCy^.-pt.m^.m^.mg), (18) 
0=1 
3 

V™ D ii = S V V* V l33 Di 3 ^ ' ~ Pt ' ~ V " m Up ' TO ^ ' m ~9 ' m U (i ) ' ( 19 ) 



0=1 
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where V is squark mixing matrix and B, C and D are loop functions defined in [27]. 
Then T^ v is given by 

r ^ = rr + rr + r ^ i + r ^ 2 

with T£" = T^\ (kl ^ k2 ^^ u ) and T B v ox2 = Y^ oxl \ {kl ^ k2 ^ v ^ u ). and T^ oxl in Eq.(20) 
the form of Eq.(4) and their corresponding coefficients Cj are given as 

c t ,i = r^(m t Vi3(C 2 d 2 - C d 3 ) - m g V 23 Cf 2 ), 
t — mi 

-4 

c t ,4 = j n(m t V 13 (C^ - C 2 e 2 ) + m g V 23 Cl 2 ) + c M , 

Ct > 5 = 2 M ^ B o + B i) - rn s m t V 23 BZ), 

mt {t-mi) 



-m g m t V 23 B\ - 



m?(t-m?) 9 (i-m})(i-miy 
x {iV 13 {B c + Bl) - m g m t V 23 B c ) 
2 ~ d 2 - 

V\ 3 C 2i + J 2"^ 13( ^24) 



t — mi t — mi 



C t ,6 = -2Vi3(C 2 e 3-Cf 2 )+ Ct ,5, 

Ct < 8 = a,/ 2 A m 2 t VM + Bl) - m s m t V 23 BZ) 
2 ~ 2 



+ ^71 ^mgmtV^BQ -r — } » 2T 

mi{t-mi) (t - mf)(t - m 2 c ) 

x (tT 13 (B C + Bi) - m a m t V 23 B c ) 

(tV 13 (Cf 2 + C d 3 ) + m 2 t V 13 (C d 2 - C d 3 ) 



t — mi 



+2V 13 C d 4 - m g m t V 23 C d 12 ) - ^—V^C^ 



t-m z t 



ctAO = 2li mt 2 A m 2 t V 13 (B^ + Bl) - m g m t V 23 BS) 
mf(t-mi) 

(iVi 3 I IK, + Bl) - m g m t V 23 B c Q ) 



(t - m 2 ) (t - mi) 

c *-i6 = r^( m *^3(c 2 e 3 - c 2 e 2 ) + m g V 23 C e 12 ), 

t-mi 

Q.18 = r^—Arn t V 13 (C d 22 - C d 2Z ) - m g V 23 Cf 2 ), 
t-mi 

Ct, 19 = — Q,20 = — Ct,5, 

csoxi.i = 4(m t Vi3 - m g V 23 )(D 22 - D 24 ) + Am t Vi 3 {D 32 - D 3e ), 
CBoxi,2 = 4(m t Vi 3 - m g V 23 )(D 23 - D 25 ) - Am t Vi 3 {D 3W - -D39), 
CBoxi, 3 = -4(m t Vi3 - m g V 23 )(D 25 - D 2e ) - 4m t Vi 3 (D 3W - D 37 ), 
CBoxiA = -4(m t Fi 3 - m g V 23 )(D 24 - D 2e ) - 4m t V 13 (D 3e - L> 37 ), 
CBoxi,5 = -4Vi 3 (L) 3 ii - D 312 ), 
CBoxis = -4y 13 (D 3n - D 313 ), 

CBoxl.l — 4Vi 3 £'3i3, 
CBoxl,8 = 4Vi3-D3i 2 , 

cbox1,9 = -4(m t Vi3 - m g V 23 )D 27 + 4m t V 13 D 312 , 
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CBoxi.ii = -4V r i 3 (£>22 - D 2A - D 34 + Daa), (39) 

CBo X i.i2 = -m 3 {D 23 - D 25 - D 35 + Ass), (40) 

CBoxl.13 = 4Vl 3 (£»25 - #26 - A»10 + A»b), (41) 

CBoxlM = 4Vi 3 (£»24 - £>26 - A»io + D 34 ), (42) 

csoxi.19 = 4Vi 3 (D 2 7 - -Dsn)- (43) 
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